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ABSTRACT 


In this work a perturbative linear response analysis is performed for the time evolution of the 
“quasi-conserved” charge of a scalar field. One can find two regimes, one follows exponential 
damping, where the damping rate is shown to come from quantum Boltzmann equations. The 
other regime (coming from multiparticle cuts and products of them) decays as power law. The 
most important, non-oscillating contribution in our model comes from a 4-particle intermediate 
state and decays as 1/f 3 . These results may have relevance for instance in the context of lepton 
number violation in the Early Universe. 
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1 Introduction 


Recently, in a series of papers [ill], Joichi, Matsnmoto and Yoshimura analyzed the effects 
of the quantum-improved Boltzmann equations. The authors developed a model of unstable 
particles and found that their Boltzmann equations generate a non-exponentially-suppressed 
equilibrium distribution function for heavy particles at low temperatures. This deviation from 
the standard picture was due to cut contributions far from the quasi-particle pole. These works 
were later criticized by several authors as a. their main point being that the definition 
of particle number is delicate for unstable particles. Redefining particle number may change 
the coefficient of the non-exponential term ||, or completely cancel it [2[. In [|3J] the authors 
argued that the naive definition needs renormalization in perturbation theory and proposed a 
non-perturbative definition for the number density. 

Non-quasi-particle properties of the quantum fields may change also the time evolution of 
the physical observables, as was also emphasized in |1], The same phenomenon was found also 
in the time evolution of bosonic condensates || |6|. In these latter cases power-law decay was 
shown up even in the linear response theory. 

In this paper we use linear response perturbation theory to investigate the long-time behav¬ 
ior of composite operators. As a specific example we consider the “quasi-conserved” charge of a 
scalar held. In this way we cannot address the equilibrium state; still, it is a physical scheme at 
which practically all systems must arrive after long time evolution. The numerical advantage 
is that all expectation values are defined in equilibrium, so that they can, in principle, be cal¬ 
culated without additional approximation. It is very important that the different contributions 
are additive, which makes it possible to identify different physical effects without entangling 
them. 

In fact there is a close analogy in the formalism describing the time evolution of the held 
condensates and the scalar charge operators; the details will be explained in the paper. In 
linear response theory the held condensates are proportional to the retarded Green’s function; 
in the case of the charge Q we have instead QF retarded Green’s function (with arbitrary F). 
The discontinuity of the retarded propagator for free case is a delta function concentrated on 
the mass shell. In the free theory Q is conserved, thus the discontinuity of QF retarded Green’s 
function is ~ 8(k 0 ) - this is the “mass shell” of the QF spectral function. In the interacting 
case, at small coupling, the mass shell still represents quasi-particles. In a similar way the 
charge, although not conserved, characterizes the system well. In analogy to quasi-particles we 
call it quasi-conserved charge. 

The shift of the mass shell is not a perturbative effect. Indeed, when we calculate the prop¬ 
agator in perturbation theory we fold IR divergences on the mass shell. One uses Schwinger- 
Dyson equations to resum these divergences, and a Breit-Wigner approximation to describe 
the shift of the mass shell. The same phenomenon appears also in the QF propagator: the 
shift of the “mass shell” causes IR divergences in perturbation theory at ko = 0. In coordi¬ 
nate space these divergences are manifested as secular terms. These, as was shown by |7|], can 
be resummed to quantum Boltzmann equations. In momentum space the same phenomenon 
shows up as pinch singularities 0B00. In order to resum them we have to consider ladder 
diagrams n, m mnum, which we will show to result in Kadanoff-Baym equations. These 
are, at lowest order [T7|], the usual Boltzmann equations. Thus the analogy of the Schwinger- 
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Dyson resummation in Breit-Wigner approximation is the use of Boltzmann equations in the 
composite operator case. 

The analogy, however, extends even more. In the propagator we find, apart from the 
broadened mass shell, other structures as well. Most important are the analytic defects (e.g. 
multiparticle cuts), which play an important role in the long-time evolution, as they give power- 
law decay mi- These are not related to the broadening, and they are perturbative in the 
massive case. We find similar defects in the QF propagator, giving power-law decay at large 
times, which are independent of the Boltzmann equations. 

This paper is structured as follows. In Section |2] we work out how initial conditions can 
be taken into account in linear response theory and what kind of long-time behavior can be 
expected. It has two ingredients: a pole contribution, which yields exponential damping, and 
a power-law correction coming from other analytic defects (such as cuts). In Section [3] the first 
contribution is considered and it is shown that this comes from the quantum Boltzmann equa¬ 
tions (Kadanoff-Baym equations). In Section [4] the off-shell analytic defects are investigated. 
Multiparticle cuts give a power-law, oscillating time dependence. It is shown, however, that 
they can be combined even at the linear response level giving non-oscillating time dependence. 
In Section |5] we summarize our findings. 

2 Long-time behavior in linear response theory 

2.1 The model 

Inspired by the discussion about the particle number, we choose a well-measurable observable: 
the charge of a bosonic field (cf. J|). The Lagrangian 

C = — lb(<9 2 + m%)\V — yt( d 2 + m 2 x )x — ^(d 2 + m|)$ — h$(^x + X Tv !) (1) 

contains a heavy charged T particle (with mass mq,) and a light charged y particle (with mass 
m x ) coupled through a Yukawa term with a scalar (with mass m$). 

Without the Yukawa term, the U(l) phase symmetry would be perfect, and the currents 

J* = d T and J* = yb d y (2) 

would be conserved quantities. For h ^ 0 they are not conserved 

dJ* =ih${ y t T-T t y), dJ x = zh$(T t y — y^T), (3) 

only their sum is. Our final goal is to describe the time evolution of the charge of T and y 

Q* {*) = (' Q ® 0)) = Tr pQ * (t) , Q* (■ t ) = J d 3 x J* (x, t) , (4) 

and similarly for y, where p is the initial density matrix. In order to avoid complications with 
the chemical potential, we assume that the total charge is zero. 
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2.2 Linear response theory 

Tracing the time evolution, starting from a non-equilibrium density matrix, is usually a hard 


task [0 . Here we follow a way where we construct the initial state from equilibrium by acting 
on the system with a time-dependent external force - after all this is the way one acts in a real 
experiment. At t = — oo we start from equilibrium, then we prepare the initial conditions by 
changing H —> H + H ext (t ) for times — oo < t < 0, finally we turn off the external force and 
start to measure for 0 < t < oo. 

This picture can be formulated in path integral representation, since at t = 
equilibrium ensemble. For an observable O(x) 


-oo we have 


(0(x)) = z 


O(x) e 


iS-i f C F 

: Jc * _ 


(5) 


KMS 


where J KMS denotes path integration where the paths are subject to KMS conditions along the 
applied Keldysh contour c. The external force is F and its explicitly time-dependent amplitude 
is denoted by ((t) (here t is the real time, not the contour time variable). 

By the end of a realistic time evolution, we deviate only slightly from equilibrium; therefore 
we expect that the corresponding state can be produced using a weak force. That means that 
in this case we can use linear response theory. If in equilibrium (0) eq = 0 then we have 


(o(x)) = -i J tiVe(f - 1') ([o(i),f(4,((i'). 


( 6 ) 


Here the expectation value of the commutator has to be taken in equilibrium, and therefore it 
is a function of x — x'. We denote 


< F) M = 0(0 <[O(O,c(0)]) e5 , e ' OF \x) = ({o(x),F( o)]>„ ; 

note that in Fourier space g^ OF \k) — i Disc Gf? F \k). 

ko 

For our concrete observable, 

Q(ko) = (Jo(k 0 , k = 0)) = g { R JoF \k)((k)\ k=0 . 


(7) 


( 8 ) 


This is very similar, in spirit, to the usual Kubo formulae [ 191, the only difference being that 


now the external force is not fixed. In principle we can imagine a lot of possible choices (e.g. a 
plausible one F ~ Hj where Hi is the interaction Hamiltonian); what is only important is that 
it should create a net Q(0), i.e. [Q,F] ^ 0. 

Having the formal result in hand, there is still a question of how to characterize the process 
and how to compute the specific numbers. In general we expect that, close to t — 0 when the 
external force was switched off, the time evolution contains a lot of transients. This is reflected 
also in the fact that Q(t) depends strongly on the choice of F. Only after a long time do we 
expect that the time dependence is independent of the initial conditions. Concentrating on this 
region we now examine the long-time behavior of an expression like (0). 
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2.3 Long-time behavior: the general case 

This analysis follows closely the one of Refs. [|, ||], including also the ( dependence. We examine 
the long-time behavior of 

( 9 ) 

where Qr is a general retarded Green’s function, ( is the amplitude of the external force. We 
assume that the analytic defects of ((k 0 ) are on the upper half-plane and the ones of Gr are 
along the real axis. Then, taking into account that using retarded Green’s functions means 
that we have to use a contour that runs slightly above the real axis in k 0 , we can extend 
the integration region for t > 0 as is shown in Fig. 0a (the semi-circle at infinity gives zero 
contribution). Under our assumptions the closed circle on the lower half-plane gives zero, while 

A 


a) 

Figure 1: Transformation of the integration contours. 



the rest picks up the discontinuity of Gr 

m = j^:e- ik °‘e(k„K(k 0 ). ( 10 ) 

The spectral function g may contain poles or discontinuities in derivatives of some order (ie. 
thresholds), but these defects form a discrete set. In between, the spectral function is analytic. 
We transform the piecewise analytic parts further. For a single interval running from —oo —a a: 

/ t#„)C(4) = e-“ j ^e-«# 1+ a)((t 0 + a). (11) 


Now we add the zero contribution of a quarter circle at infinity and a contour running up 
and down on the negative imaginary axis as shown by Fig. 0b- The closed loop picks up the 
contributions of the poles inside; in the rest we change variables ik 0 —» y and write 


Res 

c^Gloop - 


g{u) C(uj)e 


—iut 


+ e 


—iat 


7T - e yt Q{a - iy)((a - iy). 
Ztt 


( 12 ) 


The poles give exponential damping, while the rest, in general, can be an arbitrary function. 

So far we were exact, now we make an approximation valid for long times: we power-expand 
g and ( around y = 0. Since for ( the threshold a is not a special point we expect that it starts 
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with a constant. On the other hand g usually starts with Ry a , where R is a constant and 
a/0. The omitted terms are at least one power of y larger than the leading ones. Therefore 


12n e ytg ( a ~ iy ^( a ~ iy "> /^ e ytRya C( a )( 1 + £%)) = 


dy 


R({a)T(a + 1) 
2nt a+1 


( 1 + °<7 >)■ 


(13) 

The most important term, therefore, comes from the leading power behavior at the threshold. 
Finally we find for the long-time behavior: 

t^oc v- T. f t 3 ,/ 3 -iujt 1 , V- / , sR({a)T(a+ 1) e~ iat 


f(t) ™ 5Z Res lut 

uj(E poles 


+ J2 (=*=)- 

aGthresholds 


2tt 


t a+1 ’ 


(14) 


where the + sign is valid for the starting and the — is for the ending of an interval. 

An important feature of this formula is that the different terms have independent, initial- 
condition-dependent weights (values of ( at different points are independent). Constraints can 
be formulated in the form of sum rules, as for example 

Q(t = 0) = J^g(k 0 )C(k 0 ), (15) 

but it gives information on the integral properties, and not for single values which appear in 
(0- So, practically, the relative weights of the different parts in the long-time expression (0 
are arbitrary. 

This behaviour is, in fact, well known if we have a finite number of poles only, 


JV-l 

Ht) = E Cje-™ 1 , (16) 

i =0 

which can be the result of an TVth order ordinary linear differential equation with constant 
coefficients. The smallest Im uy will be dominant for long times, the other frequencies represent 
a transient behaviour. The coefficients q are free parameters coming from the initial conditions. 
For a finite number of poles, sum rules of the form 

f {n \t = 0) = £ c^-icu^, n = 0 ... N - 1 (17) 

(the superscript (n) stands for the nth derivative) can fully determine the coefficients - ie. the 
knowledge of N derivatives at the initial time is needed. In case of continuous spectral functions, 
however, detailed knowledge of the history is necessary for their complete determination. 

We can also find the analogue of the power law behaviour for the finite case. If we consider 
two oscillators working at nearby frequencies oj ± A, then 


cos(u; + A )t + cos(u; — A )t = 2 cos ut cos At, 


(18) 


the result is an average frequency oscillation modulated by a low frequency oscillation. In this 
way high frequency modes can have IR behaviour. When the number of the oscillators grows 
and at the same time the frequency difference vanishes we have 


i TV—1 

1 ^i(u+nc/N)t _ 


N 


n =0 


i- — e i<ct<7v 


ct 


_/V 1 — e ict / N 


( 19 ) 


i.e. we arrive indeed at a power-law decay. For N —> oo this will be valid for all times. 
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2.4 Long-time behaviour for quasi-conserved quantities 


The structure of g is special for an observable that is conserved in some limit, and its breaking is 
weak. In our specific case in the h —> 0 limit Q is conserved, so we expect that for small Yukawa 
couplings the following will be true. If there were no breaking at all then Q = i[H,Q\ = 0, 
therefore 

g(t, k = 0 ) = {[Q(t), F(0)]) eq = ([<3(0), F(0)]} eg = p(0, k = 0). (20) 

Its Fourier transform therefore is g(ko) ~ S(ko), as can be seen in Fig. |]a. If Q is not conserved, 



Figure 2: Schematic plot of the spectral function for a.) conserved quantity and b.) weak 
breaking. 


we expect, in general, that we have non-zero g also for other values of ko. Close to the conserved 
limit, we shall recover something like Fig ^1). The delta peak at zero k 0 will be broadened, but 
it still represents a pole near the origin. Moreover we expect to find some other structures 
(analytic defects) as well. From this spectral function, according to the previous subsection, 
we have the time dependence for long times 


0W = +0.5=^ + ..., (2D 

where T is the shift of the k = 0 pole, is the position of the threshold with the smallest a. 
The coefficients Q i, Q 2 and 0 depend on the initial conditions. 

One could think that power-law terms dominate for really long times. This is true; still, 
there are effects which make them less relevant. One is that usually O ~ m (a mass scale in the 
system), and so they are fast oscillating functions. Thus they can hardly influence physics on 
scales larger than m -1 . The O ~ 0 case is exceptional. This, however, shows up at higher order 
in perturbation theory (or comes from non-linear effects - see later) and Q 2 is proportional to 
some powers of the coupling constant h. Thus Q 1 ^ Q 2 (unless we have very special initial 
conditions), and the first term will be larger for a large period of time. Analyzing the functions 
we find that the two terms are equal at time t, where 


t <Tt < 1.582t, 


r = In 


Q 1 / ol + 1 


a-\-l~ 


( 22 ) 


Since Q1/Q2 and T are polynomials of h, the first term is the leading one up to times Ft ~ 
—Clnh, where C ~ 5-10 typically (for example if Q2IQ1 ~ h 2 , a + 1 = 3 and T ~ h 2 then 
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C — 8). During this time the charge drops a factor of e~ rt ~ h c , i.e. for small couplings the 
power law regime may be unobservable. A hint for the presence of such a regime, however, can 
be seen in numerical simulations 
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On the other hand, if h is not small, then the exponential behaviour may be suppressed. 
In fact, only in quasi-particle picture, when T <C other mass scales, do we expect exponential 
damping. 


3 The damping rate 


In the previous section we outlined the general long-time behaviour for a quasi-conserved quan¬ 
tity such as the particle charge. The two relevant effects found are broadening of the delta 
peak at k 0 = 0, and analytic defects at k 0 = 0. In this section we make some remarks on the 
broadening effect. 

In our specific example at h = 0 the charge is conserved, and we assume that for small h 
it still represents a pole on the second Riemann sheet close to the origin. For small k o we thus 
expect that the time-ordered or retarded propagator (at least the pole contribution to that) 
can be approximated as 

G( QF \k 0 ) = (TQF)(k 0 )^—^—, for k 0 - 0, Im k 0 < 0. (23) 

x ' k 0 + zL 

This form suggests that the broadening effect must show up as IR divergences in the pertur¬ 
bation theory. Here, namely, we can calculate only h -power corrections to the leading result, 
therefore the appearance of T in the denominator means that in perturbation theory we find 


G (QF \k 0 -h. 0) 


G ^ f-iT 

fc o„= 0 


kn 


(24) 


The signal of this behaviour is the ~ l/fog divergence of the perturbative expression for small 
(but not zero!) k 0 . 

Working in coordinate space the k^ n divergences turn into t n ~ l time dependence. These are 
the secular terms that appear in all orders in perturbation theory, and that can be resummed 
according to [|7|] as the quantum Boltzmann equation. 

Using perturbative methods in Fourier space, the origin of these IR divergences is the 
pinch singularities [|, |9jj. To see this consider a diagram contributing to the QF retarded 
Green’s function with a 2-particle intermediate state, as can be seen in Fig. [|. The momenta 
of the propagators G\ and G 2 , because of energy-momentum conservation, are q + k 0 and q, 
respectively. If these propagators are retarded and advanced with the same mass, then for 
ko = 0 the poles of the two propagators approach the real axis from above and from below, 
thus “pinching” it, giving the divergent result GR(q)GA^q) ~ {d(q 2 — nr 2 )) 2 ||, |J. For finite k 0 
we do not expect divergences, but the ko —■> 0 limit must be infinite. In Appendix |A] we indeed 
find that, at the k 0 —■► 0 limit, we have a singular behaviour for the products 


G r 



G a q~ 


ig(q) 


d (q) 

2q 0 k 0 ' 


Gr/a 




d (q) 

2q 0 k 0 


G rr (q), 


(25) 
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Figure 3: IR-divergent (ladder) diagram. The shaded insertions can represent rungs or self¬ 
energy insertions. 


and the inclusion of G#(g + k/ 2) yields an extra d(g)/2go&o factor, while the inclusion of 
additional GA(g — k/2) yields —d(q)/2q 0 k 0 . Therefore a ladder containing n rungs or self¬ 
energy insertions may yield a k q h IR divergence, corresponding exactly to our expectations. 

We can also argue that more particle intermediate states do not give pinch singularities 
in the massive case. Consider a 3-particle intermediate state with momenta p, q and — p — q 
{k = 0 case). We expect an IR enhancement if the sum of two momenta is almost zero; let us 
take q ~ p + q, i.e. p ~ 0. Then the propagator with momentum p gives 1/m 2 , where its mass 
is m, the other two propagators are replaced using Appendix |A]. One finds, concentrating on 
this region, 

^4 ^4 Q(q) 1 for\ 

dpdq—pr^—. -—-—?■ ( 26 ) 


m- 2p 0 q 0 - u/ +q + u, 


This expression is IR-safe, as can be seen for example by integrating over the IR region |po|? |p| < 
e. Therefore we cannot have any 1 /k 0 enhancement from this insertion. This argumentation 
can be extended to more internal lines. 

There exists another possibility when there is an IR sub-divergence, coming from a self¬ 
energy insertion on an internal line. Since, however, internal lines are equilibrium propagators, 
these pathologies must finally cancel [§]. 

To resum the ladder diagrams we perform point splitting in a symmetric way: 


(M x )) = - id/)- (T(x)T f (|/) + T t (r/)T(a;)) = lim(id£ - id/)iGZ(x,y), (27) 


y-*x' 


and similarly for y, where we used R/A formalism to define the propagator [23]. The Q 
propagator corresponds to the boxed part in the diagram of Fig. |3j. In Fourier space 


(M k )) = J 77 ^h ( 2( ? + k )»iGrr(q + k,q ). 


Introducing Wigner transformation 


W(g,X) = jd l ue-^iQ (V + |,V-|), W(<f, *) = iQ(g + q - |) 


we can rewrite it as 


= 


(Pq 

W)' 




(28) 


(29) 


( 30 ) 

















A ladder diagram is a recursive structure, thus there must be a kernel which relates a 
ladder containing N + 1 and N insertions. Keeping only the pinch singular terms we can write 
iQ N+} (p,p') = 1/2 k 0 q 0 (K *iQ N )(p,p'), where q — (p + p')/2 and k = p — p' (cf. (^)). The sum 
for all ladders Q = J2nG n then satisfies iO{p,p') = l/2k 0 q 0 (K *iG)(p,p'). With the Wigner 
transformed functions fl29P , this can be written as W(q, k) = l/2koqo(K * W)(g, k), or, in real 
space (performing a Fourier transformation with respect to k), 


2 iq 0 dtW{q,t) = {K*W){q,t), 


(31) 


a hrst order differential equation. 

This differential equation is unique up to leading order in the k expansion. Let us assume, 
namely, that there is another equation 

2iq Q dtW(q,t) = (. K’*W)(q,t ) (32) 

satished by W. Then inserting the solution of ( |?I[ ) into this equation we hnd 

(K - K') * W = 0. (33) 


This is an identity that must be fulfilled for K and K'. This shows that K = K' on the space 
of the solutions of (pl|). 

Using this property we can use differential equations of the form of (|3lJ) with different 
origin in order to resum pinch singularities. And exactly of this form are the Kadanoff-Baym 
equations |2(| in linear response. For the free theory, namely, the propagators satisfy, because 
of the Schwinger-Dyson equations, 


(dl - d$)g(x, y) = 0 =*► 2 iq 0 d t W(q, t ) = 0. (34) 


In the case of interactions the right-hand side depends on Q (or W), but in the linear response 
theory it is linearized and has the form K * W with some kernel. We therefore recover (|32l). 
The uniqueness of this equation then implies that the Kadanoff-Baym equations represent an 
adequate tool to resum the pinch singularities of ladder diagrams in linear response theory. 

So finally both the coordinate space method [|7J and the momentum space method outlined 
above support the statement that the resummation of l/kfi IR singularities of the perturba¬ 
tion theory yields quantum Boltzmann equations (Kadanoff-Baym equations). Therefore T 
appearing in (f23[ ) must be calculated using Boltzmann equations. 

In our model to lowest order we have to take into account the decay processes, which have 
a matrix element h. Therefore the Boltzmann equations are 


2 qd x N*{q,X) = h 2 


■ dH 

Wy 


*(q-£) (-/V w + l)jV x IV® — Aw(l + N x )(l + N®) , (35) 


and similar equations for x and < F. In the integrand we denoted N^(q,X) —> , N X (£,X) — > 
N x and N^^q — £,X) —> N®. The possible values of q 0 are q 0 = ±tu q (the dispersion relation). 

We solve the equations in the relaxation time approximation. We rewrite (|3(J) using our 
condition that there is no net charge in equilibrium 


Q 9 (t) = 2qoQ*(q)SN(q,t), 


(36) 
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where SN is the deviation from equilibrium. For SN the linearized equation reads 


2 qd x 6N* 



y 

9 X 0 


'8N*(l + N X + N*) + SN X (N* - N + 5N*(N* 



( 37 ) 

We multiply this equation by g*(q) and integrate over q. We can realize that, because is a 
real field, 


9%{x, y ) = g%{y, x ) =► Gf r (p, p) = Q%{-p', ~P ) 


8 W*(q,k) = SW*(-q,k). (38) 


Since SW = g^SN ]TjJ, this symmetry renders the last term of 0 zero after integration. We 
introduce 

r ^>= wJw? ^ e)s% ~ i](i+m)+N[q ° - 4,) = ~ 2 y w - 

r x( 9) = ^/(01 /(«)/(<?-0(1 + m«) + JV(9o - «o)) = (39) 

where N(x) = (e^ x — 1) _1 is the Bose-Einstein distribution, and write 


dtQ 9 (t) 



-2q Q 5N*( qi t)gV(q) + 2q 0 SN x (q,t)g x T x (q) 


(40) 


The relaxation time approximation assumes that T is only weakly momentum-dependent, so 
we can factorize the integral, and we finally obtain 


d t Q*{t) = -r*Q*(t) + r*Q*(t), 

dtQ x (t) = —T X Q X (t) + r (t). ( 41 ) 


This approximation satisfies the requirement that Q' v + Q x = const., which is zero (the equi¬ 
librium value). Then the solution of the above equations is 

Q*(t) = Q*(0)e- (rll/+rx)t , Q x [t ) = -g' I, (0)e- (r " /+rx)f . (42) 

For m<a S> m x , we hnd 

7 o 7 o 

r* = ^ (l + 2 N(^)) , r x = —- (3MN(M)(1 + N(M)), (43) 

107rrrnp V 2 / lo7rm x 

where M = m%/2m x . 


4 Contributions coming from analytic defects 

The other important ingredient of the long-time evolution is the power law tail coming from 
analytic defects of the spectral function. These, in general (cf. Section Q), yield the time de¬ 
pendence Acos(ttt + (ft)/t a+1 . The value of and a depends on the theory, while the amplitude 
A and the phase <p come from the initial conditions. In this section we try to compute and 
a for some special cases. 
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4.1 Two-particle cut 

Let us first consider a diagram where there is a two-particle intermediate state, as is illustrated 
in Fig. [|a. We represent it as two momentum-dependent vertices denoted by Ti(q,k) and 



Figure 4: Two- and multiparticle intermediate state 


r 2 (g, k ), their product being denoted by Y = The intermediate particles have masses m\ 

and m 2 , and we treat it as a cut diagram (with Keldysh indices 21 minus 12). Its contribution 
then is 


C 2 (k) = 


■ d A q 

Jiiy 


r(?, k) \Gl\q)G?(k -q)- G?(q)G?(k - q) 


= 47T 2 


d 4 q Y{q,k) 


[(1 + N(q 0 ))(l + N(ko ~ go)) - N(q 0 )N(k 0 - q 0 )] 

Qo — ui 2 ) — S(ko — qo + ^ 2 )) • 


(27t) 4 4u;icu2 

x (%o — ^i) — + <^i))(^(^o 


(44) 


The analytic defects of this function are thresholds, where the phase-space, restricted by the 
constraints of the integration, vanishes. There are 4 such points in this case: ko = ±mj ± m 2 . 
We restrict ourselves for positive k 0 and denote k q 1 * = mi +m 2 and k^ = mi — m 2 (assuming 
m 1 > m 2 ). At h 1 ' 1 ' 1 a cut starts, at k^ it stops, so for the threshold behaviour we choose 
ko = k^ + y or ko = k^ —y. Since for y —> 0 the phase space vanishes (q —* 0), we can expand 
everything to second order with respect to |q|. 

In the Ah 1 ) case the relevant combination of the delta functions of (fnj) is 

S(q 0 -uJi)S(k 0 - g 0 -a; 2 ) =>■ A: 0 = a;i+a; 2 . (45) 


Expanding u>i + o; 2 around k^ = mi + m 2 to second order yields 

V = 


1 1 
2mi 2 m 2 


q 2 , 


(46) 


so that |q| vanishes as y 1 ^ 2 . Assuming that Y does not vanish in this limit (as the best case), 
and we substitute u;i )2 —> mi )2 , k 0 = mi + m 2 and q = 0 where it is possible, we find 


C ~ (1 + A”(mi) + N{m 2 )) JS (y - 


d 3 q s ( / 1 

\2mi ' 2 m 2 


+ 


q‘ 


M 2 


(47) 


which can be seen after performing the integral. This is the same as the zero-temperature result 
21]], only the amplitude is modified by (1 + N(mi) + N{m 2 )). 
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The fc< 2 ) case goes along the same lines. We choose the combination 

5{qo — u)i)5{kQ — q$-\-U 2 ) => k$ = u)\ — uj 2 => 2/=^-m- 2m ) 




and therefore 


C ~ (N(mi) - N(m 2 ))y l /\ 


( 48 ) 


(49) 


0 - 


This is a finite temperature effect, the threshold of the Landau cut 

For our concrete example described by the Lagrangian ([!]), only equal-mass two-particle 
intermediate states are allowed. In that case the first threshold is just a zero temperature one, 
starting at 2m* where m,; e {m, M,ms}- The second case, however, is not present, since the 
delta constraints then give S(ko)- So there the important cut contribution has O = 2m* and 
a = 1/2. That is we obtain an oscillating contribution which has presumably negligible effect 
on real physical processes. 


4.2 Four-particle intermediate state 

Four-particle intermediate states may contain cuts, and this will be analyzed in the next sub¬ 
section. Here we concentrate on other possible analytic defects. 

On the same footing as before we can write up the contribution coming from the four-particle 
intermediate state 

cm = ..4hL s(k - E 9 )r(k, {,}) [(1 + + Ay - Ay.Ay (so) 

where N t = N(q 0i ). We change q —> —q in the first term and write C±{k) = C±{k) — Cff (k) 
where 

cm = s(k ± Eq)m m}m.N 4 (sy 

In fact, this is a general form; we can write also for the two-particle contribution ([44]) C 2 (k ) = 
C 2 (k) — C 2 (k ), where 

c m = j ^i^ i ^±Y.q)nk,m})N 1 N 2eiQ2 . (52) 

We now make a simplification, and assume that for some reason T is factorized as T(k, {g}) = 
r ( 1 ’ 2 )(fc, q 1: g 2 )r ( 3 , 4 )(A;, q 3 , g 4 ). This means a loss in generality; but we hope that the main 
consequences still remain true. We introduce an auxiliary variable Q and write 

Ct{k) = jJ (2t!) 4 (2tF) 4 + ?2)r (1,2) (fc, —qi, -q2)NiN 2 Q l (q 1 )Q 2 (q 2 ) 

X / ( 2 n ) 4 ( 27 t ) 4 6 ( k ~Q + , i 3 + g 4 )r ( 3 ’ 4 ) (fc, -g 3 , -g 4 )^ 3 ^ 4 ^ 3 (? 3 )^ 4 (g 4 )-( 53 ) 

On the right-hand side we recover the convolution of two C 2 contributions; in real time they 
are simply multiplied. The complete four-particle contribution can be written as 

C A (t) = C+(t)C' 2 + (t) - C 2 (t)C'-(t). (54) 
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At zero temperature N(q 0 ) —> — 0(— q 0 ), thus contains only the positive energy cuts and 
C+C't oscillates with frequency Cl + Cl'. At finite temperature, however, the negative energy 
cut appears as well, and, for the same threshold values, we have 


^iQt _j_ ^ ^ g—iQt 

W 2 W 2 


A_A' + + A + A'_ 

T 3 


+ oscillating. 


(55) 


Diagrammatically it comes from the diagram where two lines are incoming, two are outgoing, 
or vice versa. Kinematically, with the external line, it is a 2 —> 3 or 3 —> 2 process. Also from 
here we see that this contribution must be absent at zero temperature. We can assess the 
temperature dependence by the typical factor for 2 incoming and 2 outgoing internal lines, 


(1 + N(rrii))(l + N(m 2 ))N(m 3 )N(m^). 


(56) 


For low temperatures the contribution is therefore Boltzmann-suppressed. 

We can also have cancellations at finite temperature, if the + and — parts in (|54p cancel 
each other. C + and C~ are related as (writing out explicitly the dependence on the T vertex 
function) 

C + (T,-k) =C~(t,k), where f (k, {q}) = T(-k, {-q}). (57) 


Cancellation occurs, if T = ±r for small k. For generic F at finite temperature we do not 
expect such a symmetry, therefore a net effect will remain. 

How can this phenomenon show up in a diagram? If we calculate Q(ko) for small k o, we 
encounter ladder diagrams, as is shown in Fig. [|. The ending of the ladder (the triangle at the 
right-hand side) contributes to the numerator (the constant C in eq. (|23|) ). Here can sit the 
four-particle intermediate state of this section, see for example Fig. |5|. The fc 0 ~ 0 behaviour 


k 0 



Figure 5: The realization of a four-particle internal state. The shaded oval denotes ladder 
resummation. 



of the charge is therefore 


Q{k o —> 0) ~ 


c 2 + #c 4 (fc 0 )h 4 c 2 + #c 4 (fc 0 )h 2 c 4 

k 0 + *r 


(58) 


The pole and the k 0 = 0 defects are therefore separated and do not affect the behaviour of each 
other. When computing the numerical coefficient of the k 0 = 0 defect we substitute k 0 = 0 in 
the numerator, which means an enhancement factor of 1 /T. 

Four-particle initial conditions can also come from the quadratic response (as a product 
of two two-particle operators). Having the same effect, it is natural to assume that their 
contribution is also in the same order of magnitude, i.e. <(4 ~ Cl- Since C •C 1 (linear response 
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regime) we find suppression in this terrrif]. To have only one small parameter in the system we 
will assume that (2 h 2 . 

Summarizing, the 1/t 3 contribution of the four-particle intermediate state is suppressed by 
several effects compared to the Boltzmann pole, ft has a perturbative suppression of h 2 ; it is a 
finite-temperature effect, so that it is Boltzmann-suppressed at small temperatures; it has two 
contributions with opposite sign; and finally it is a multiloop effect, which probably means also 
a numerical suppression. 

4.3 Multiparticle cuts and other defects 

The determination of the threshold behaviour described in the case of two-particle intermediate 
states can be generalized to multiparticle intermediate states (cf. Fig. f|&). We expect thresholds 
when the external k 0 plus the N internal particles represent a decay process (1 —> N or N —> 1), 
since scattering processes do not have cuts. There always exists a possibility, when all of the 
internal particles are outgoing, so that = J2i mi is the threshold value; this is also the 
zero-temperature case. At finite temperature an additional threshold is possible (for positive 
ho), when one particle is outgoing, the rest are incoming. Then we obtain a threshold at 
k q 2> = m 1 — Y^i= 2 rn ii if it is positive. 

Near the threshold, in the massive case, we can perform the same expansion as in ( |46|) , and 
also there is a momentum-conserving delta function S(J2 q?;)- Therefore what remains is 

/ d 3 q,. ..d 3 q m 5 ( 9 - £ |h) S ~ y ^-\ (59) 

The possible values for thus are trk and m\ — (if it is positive), in both cases 

a = 3(Af ~ 1) — 1 (which is also the zero-temperature result). 

Similarly to the two times two-particle cut case, multiparticle cuts can also be multiplied 
in real space to generate new defects. This means, however, a faster decay than 1/t 3 . 

5 Conclusions 

In this paper we tried to describe the time evolution of the charge operator Q in a complex 
scalar system, starting from an initial state not too far from thermal equilibrium. The initial 
state was prepared using an external action that modifies the time evolution for a given period. 
We applied linear response theory in the external action to describe time evolution. As a 
computing method we used real time perturbation theory. 

The long-time behaviour of Q is determined by the QF spectral function in the linear 
response theory, where F is an arbitrary operator ([Q,F] 7 ^ 0). If Q were conserved, there 
would be g( QF) (k 0 ) ~ 5(ho). For the non-conserving case two types of contributions can be 
identified (cf. IS! ). One is coming from the broadening of the S(ko) pole of the conserved 
Q-case, yielding exponential damping for long times. In perturbation theory it manifests itself 

2 This hierarchy, however, does not necessarily hold: wild initial conditions can violate it, then the 1/t 3 time 

dependence also has a coefficient of 0(1). 
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as pinch singularities in the ladder diagrams, giving 1/fcg type of IR divergences at any order 
in the coupling. Resummation of these singularities results in linearized Kadanoff-Baym equa¬ 
tions. Working in real space, 1 //eg singularities correspond to the secular terms t n_1 , and their 
resummation also leads to quantum Boltzmann equations [j7j. 

Analytic defects in g^ F \ko) also contribute to long-time physics: it gives an / t a+1 

oscillating power-law decay, where is the position of the defect, a is the index of g near to 
the defect. For multiparticle cuts we have = J2i rn i or m* — X)j>i m i (if it is positive: a 
multiparticle Landau damping), and index a = 3(At — l)/2 — 1, where N is the number of 
intermediate particles and m; are their masses. These contributions are oscillating with the 
particle mass, so that they can hardly influence physics on IR scales. At finite temperature, 
however, these cuts can be combined, and from their superposition the oscillating behaviour 
can drop out. The lowest order diagram giving this result contains 4 intermediate particles and 
yields 1/t 3 type of decay. 

So finally, starting from a generic initial state characterized by F, we expect that the time 
evolution of the charge for large times follows 

Q(t) = Qie~ rt + ^ + oscillating. (60) 

t 6 


The relative weights of the Boltzmann pole and the cut contributions contain free parameters, 
depending on the initial state, and they contain perturbative factors. The pole is present 
already at the free level, we thus expect Q\ = 0(1) + perturbative corrections. The amplitude 
of Q 2 is reduced by several effects; for example it is suppressed by h 2 and by the Boltzmann 
factor for small temperatures. Thus, for a generic situation we expect that exponential damping 
dominates the time evolution for intermediate times and only after a long time will the power 
law take over (for a numerical simulation example see, (22)]). Still, because of its slow decay, it 
can be relevant in the explanation of small deviations from equilibrium (such as baryogenesis). 
This should be examined in the future. 

Another task for the future is to include non-linear effects and see how the full Boltzmann 
equation gets modified due to the presence of cut contributions. 
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A Appendix: pinch singularities 

The pole structure of the retarded and advanced propagators is represented by 
Dr / a ^ = (g 0 ± ieY - oA = 2ZJ 


g 0 — u ± ie g 0 + u ± ie 


(61) 
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where u;(q) is the dispersion relation. They have poles at ±u; — ie for the retarded Green’s 
function and ±u; + ie for the advanced one. Their product thus has poles at ±u; ± ie , which 
pinch the real axis for e —> 0, resulting in double poles which cannot be regularized. We do not 
expect such problems for products like DrDr or DaDa- 

In our case the momenta of the retarded and advanced propagators are offset by k. For 
finite k we do not expect singularities, but the k —> 0 limit must be divergent. 

First we examine the D R (q + k)D A (q ) product for small k = (. k 0 , 0 ). We will use the identity 


1 ( l 

~ = V - 

x + ie \x 


md{x) 


where V means principal value, so that 


Dr/a(q) — 


For the product we obtain 


1 

2 a; 


V 


Qo ~w t 


=F i^d(q 0 — to) — {a; —► —a;} 




i7rV tk) 

+ir 2 5(k 0 )5(q 0 — uj) + {a; —> -a;} 


(62) 


(63) 


+ non-singular, (64) 


the non-singular terms having finite limit at k —> 0. Rewriting this expression we obtain 


. k . k. ( 1 

DR[q+ 2 )DAq -2 )= ' P {2^ 


Q(q)+ °(£)' 


(65) 


where g is the free spectral function; o(l/x) means that lim aw0 ( a; o(l/x)) = 0, which is trivial 
for non-singular terms, and true for the delta term since xS(x) = 0 . 

In a real model we can write for the propagator G(q) = d(q)D(q), where d is the Klcin- 
Gordon divisoij] ||. 

We can obtain (^) in a different way, using i(Gn — Ga) = Q and that GrGr is non-singular. 
We write 

Gr(q + -j)G A (q — —) = G R {q + —) ig(q — —) — G R (q + k)G R (q). ( 66 ) 

The second term is finite for k —> 0, while in the first term at k = 0 we put a propagator on 
mass shell, which is clearly divergent. The numerator of G^(g + k/ 2) for small k reads 

k 2 

ql + Qoko + ~ w q + k /2 —i " 2go^o + ^q_ k / 2 — w q +k / 2 —>” 2g 0 ^o — 2k • du ;, (67) 

where we have power-expanded o; q ± k / 2 . Therefore 


/ k / k d(q) 

Gr(i + 2>Ga(9 - 5 ) - 2qoka _ 2k . doJ 


ig(q) + 


3 For fields with multiple mass shell we obtain a sum of these terms. 


( 68 ) 
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where d(q ) is the Klcin-Gordon divisor. It is also clear that any other Gr factor yields the 
same d(q) / (2q 0 k 0 — 2k • du) contribution. Additional Ga factors yield — d(q)/(2q 0 k 0 — 2k ■ du) 
contribution, as can be seen by using the above formulae for Gr = —ig + Ga- 

In the paper we use k = (ko, 0) momentum, but all the formulae can easily be generalized 
to finite k simply by 2q 0 ko —> 2q 0 ko — 2k • du. 

Pinch singularities occur in other products as well. In the R/A formalism [23] the propa¬ 
gators are Gr , Ga and G rr = (G*i 2 + G 2 1 )/2. Since G rr ~ g also puts the propagation on the 
mass shell, we find 

Gr/a(q + ^ )G rr (q ~ = 2 %kj~* rr ^' (^) 


On the other hand G rr (q + k) ■G rr (q ) ~ 5(Ao)> therefore it yields an o(l/k 0 ) contribution. So in 
the two-particle intermediate states the possible pinch enhancements are coming from GrGa 
and G rr GR/A only. 
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